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Abstract. We address the problem of two interacting atoms of different species 
inside a cavity and find the explicit solutions of the corresponding eigenvalues and 
eigenfunctions using a new invariant. This model encompasses various commonly used 
models. By way of example we obtain closed expressions for concurrence and purity 
as a function of time for the case where the cavity is prepared in a number state. 
We discuss the behaviour of these quantities and and their relative behaviour in the 
concurrence-purity plane. 
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The system of two two-level atoms (TLA) inside a cavity has attracted considerable 
attention, both because it has become experimentally feasible and because it is the 
paradigm to study the evolution of entanglement under decoherence. This combination 
is remarkable, because entanglement is a central resource and decoherence the major 
impediment for quantum information processing |T]. The relation between concurrence 
and purity of the central system yields the simplest access to the problem. 

Different models of two identical TLA as a central system coupled to a cavity mode 
in resonance with the atomic transition as environment have been studied (2j |3] HI [5] . 
In this paper, we show that one can define a wider class of such systems that remains 
solvable in closed form and includes the above mentioned cases. Specifically we consider 
atoms with different coupling to the cavity mode, different detuning and include dipole- 
dipole as well as Ising interactions between the atoms. 

We show that the total number of excitations is a conserved quantity. Using 
the basis in which the corresponding operator is diagonal, the Hamiltonian will be 
transformed to block diagonal form, with maximally 4x4 blocks. Interestingly we 
could use a special case of this solution to construct an exactly solvable relativistic 
model [!(]] with three degrees of freedom, namely a Dirac oscillator [TT] coupled to an 
isospin field. 

Note that models with different interacting TLA and a single excitation on a 
continuum of modes have been solved P, [7J E] using the pseudomode approach [9] which 
in those cases results in a single mode with losses. While quite similar, the loss term 
violates the conservation law we use and thus these are not particular cases of our model. 

In order to focus on a particular new aspect, namely the interplay of Ising and 
dipole-dipole interaction, we shall choose an example where other features of our model 
are simplified. Thus we shall apply the closed solution to study the time evolution of 
concurrence and purity of two interacting TLA with equal coupling and zero detuning 
but arbitrary dipole-dipole and Ising interactions. The interaction free case basically 
provides the borders of the evolution if we look at the interacting problem in a 
concurrence-purity {CP) diagram, a third boundary being provided by the relative 
strength of the two interactions. 

Particular cases of the general model, for which solutions are available, should 
be experimentally feasible in cavity QED [121 03] • While dipole-dipole interactions 
commonly appear in QED, an Ising interaction might be simulated as proposed in 
references [HI [151 HB]- Wether a particular case, such as the one we discuss, will actually 
be measured depends on specific difficulties in forming the initial state, as well as the 
amount of interest such a case may arouse. Some such cases, including initial coherent 
states, will be studied in a forthcoming paper [TT] . 

Consider the Hamiltonian for two TLA coupled to a cavity mode and set % = 1, we 
use the rotating wave approximation and work in the interaction picture so we end up 
with 

H=T, fao? + 9i (a + aV<?>) } + 2k (<&a® + a^a^) 
i=i 
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+ M 1 ^ (i) 

where 8j is the detuning of the corresponding atomic transition frequency from the 
frequency of the cavity mode which does not appear due to our choice of the interaction 
picture. 

<?j is the coupling to the mode, k and J are the strengths of the dipole-dipole 
and Ising interactions respectively. We use the standard definitions of creation and 
annihilation operators for the harmonic oscillator (a, a') and for the raising and lowering 
operators a± = (cr x ± ia y )/2, with the Pauli matrices (a x , o~ y , a z ). 

The operator I = o)a + 1/2 (a^ + a^^j provides an additional constant of motion 
and it can be interpreted as the number of excitations in the system. Clearly [H, I] = 
and in general this is the only commuting observable of this problem. Therefore we 
choose the following basis for which I is diagonal 

\<fc ) ) = \n + l)\ — ) I0^) = |n)|-+) 

l^ B) > = |n>|+-> |0i" ) ) = |n-l)|++). (2) 

Here \n) describes a state of n photons in the cavity, |— ) and |+) describe the ground 
and excited states of a TLA respectively. For any given n they satisfy the relation 
I\(f>j^) = n|0j™' > ). In this basis H is a block-diagonal matrix and each block is a 
4x4 matrix with elements = H^K Explicitly, one has 

/ J — 5i — 5 2 92\fn + 1 giVn + 1 \ 

H (n) = 



(3) 



g 2 Vn + 1 S 2 - 5 1 - J 2k giy/n 
gWn + 1 2/c 81 — 62 — J g2\fn 

V g ly /n g 2 y/n J + 8 1 + 8 2 J 

For n — 0, the basis is reduced to the three states |1)| ), |0)| — h) and |0)|H — ). For 

n — — 1 it is reduced to one single |0)| ). This single state is stationary and represents 

the situation where both atoms are in the ground state and there are no photons in the 
cavity. 

Solving the resulting eigenvalue problem implies diagonalizing each block of the 
Hamiltonian. In general the characteristic polynomial for the eigenvalues leads to a 
depressed quartic equation with eigenvalues: 



+ ^v /2i? - f/ + 7& ^ = 1 > 2 
Ef = (4) 

^+ { ^PR-U-^ J =3,4, 
where we used the following definitions: 

P = (6\ - 81 + (n + 1) (g\ - gfj) (8f - 8\ + n (g{ - g*)) 
+ J 2 ((2n + 1) (gl + gf) - 2 (6f + 8 2 2 ) + J 2 - 4/t 2 ) 
+ 2J (g 2 ^ + g 2 2 8 2 + 2(2n + l) Kgi g 2 ) 
+ Ak (8i + 6 2 ) (gig 2 + k (<Ji + 6 2 )) 
Q = 4 (gl8 2 + g 2 2 8 x + 4J (k 2 - 8^) - 2(2n + l)K 9l g 2 ) 
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R = 2/3 ((2n + 1) (g\ + g 2 2 ) + 2 (<5 2 + 5% + J 2 ) + 4k 

n 2 - R 3 4P /? 2 

S=2PR + ^L T= T + T 

17 = (S + VS 2 -T 3 ) 1/3 + (S - VS 2 -TA 1/3 . (5) 



The eigenvectors before normalization read as 

In 



v\ n ) = (Ef ] - 5l -5 2 -J)x 

{{Ef ) + j) 2 ~n{gl + gt)-{5 l -5 2 f~^ 
- 2n ((^ + J) g\ + {5 2 + J) 5 2 + 2rc 5l £ 2 ) 

v/n+T^! (^j n) - 5i - 5 2 - J) 
+ <?2 ((if } - 5i) 2 + n (gf - gt) - (5 2 + J) 2 ) ) 
V5i+Tf2«p 2 (£j n) -6 1 -5 2 -j) 
+ Si ((if } - 5 2 f + n (gi - gf) - (5, + J) 2 ) ) 



(n) 



(n) 



«4j = V /'(/' 



+ 1) {2g 1 g 2 {Ef ) + J) + 2« (<? 2 + ^)) (6) 
and the orthogonal transformation which diagonalizes the Hamiltonian is given by 

T/(n)_..(n)// v (n)2\V2 
v j,k — v j,k I v l,k J 

By way of example we now treat a special case where we calculate the entanglement 
and purity of the pair of atoms considering the cavity mode as environment. For this 
purpose it is convenient to start with product states of cavity and central system 
functions. We restrict ourselves to a definite value of the observable I and choose a 
number state for the cavity i. e. 

|^o) = \n) (cos(a)|-+) + sin (a)|+-)) . (7) 

In the same subspace of fixed eigenvalue of I, one could also use the state \n + 1)| ) 

or | n — as initial product states. This type of initial state guarantees that the 

evolution stays confined in a four dimensional subspace. The time evolution of the state 
vector under the Hamiltonian (JTJ, can be written as 

!*(*)> =Es l w (OI# ) > (8) 
i=i 

with the following coefficients 

4 n) = £ V#>e"<« (Vff cos (a) + vff sin (a)) . (9) 
k=i 

For readability, we shall omit the time dependence in the coefficients, = B^\t). 
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Starting from the density matrix of the whole system g(t) = \ty(t))(^(t)\, we take 
a partial trace over the cavity degree of freedom to compute the reduced density matrix 
of the two TLA, given by 



( \B { r ] \ 2 




V 





\Bt ] \ 2 







Bl 



(n) 



B 



(n) 



{BP) m B<? 
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(10) 



The purity P = Tr p 2 measures the entanglement between the central system and the 
environment, i.e. the decoherence of the two TLA and we find 



P = \B{ n) \ 4 + | J Bi n) | 4 + (l " |£i n) | 2 ~ \B. 



Ml 2 



(11) 

The concurrence [18] is used to measure the entanglement between the atoms. It 
is defined as C(p) = Max {0, Ai — A 2 — A 3 — A 4 }, where A.,- are the eigenvalues of 

in non-increasing order. In our case the concurrence is given 



r (l) a (2) 

y y 



1/2 



by 



C(p) = Max{0, 2|4 ra) || J Bf } | - 2\B[ n) \\B { ; i) \} 



(12) 

Some interesting features can already be inferred by inspecting ( TTTj) and (|12l) . For 
n = we have B^ = and the purity has a minimum value of 1/2. As for the 



concurrence one can note the absence of entanglement sudden death [19, 20J in that 
particular case. 

Now we specialize in the symmetric case with equal couplings to the cavity, 
zero detunings, but allow both types of interactions between the atoms. With these 
restrictions we are able to find explicit solutions in the time domain. Using the 
definitions 



UJ,, 



0n 



'An 



An 2 + An 



In 



6n 2 + 6n + 2 
An 2 + 4n + 1 ' 



(13) 



V An 2 + An + 1 ' 

where u n is a frequency closely related with the eigenvalues of the Hamiltonian ([1]), and 
the time-dependent functions 



F(t) 
G(t) 



2n+l 



1 + sin (2a)) sin 2 (co n t) 



[k — J) cos ((J + 3n)t) sin (u n t) 



+ 



sin ((J + 3fc)t) cos (uj n t), (14) 

we find the following solutions for the purity and the concurrence as functions of time 

C(t) = Max{0, v^sin (2a) - F(t)) 2 + cos 2 (2a)G 2 (t) - f3 n F(t)} 

P(t) = l-2F(t)+ ln F 2 (t). (15) 

In figure [1] we present these solutions for the situation where there are no initial 
photons in the cavity, namely n = 0. We present three cases: in red and blue for 
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Figure 1. (Colour online) Concurrence and purity as function of time for n = and 
an initially empty cavity. The red curve corresponds two non-interacting atoms with 
an initial state determined by a — 7r/4 (See eq. ([7])) i.e. a maximally entangled pure 
state. The blue dashed curve shows the behaviour for non-interacting atoms with an 
initial pure, but not maximally entangled state with a = n/20. In black, the curve 
for two interacting atoms with the same initial state as in the blue dashed curve and 
K = 1.5 and J = 0. 

noninteracting atoms with initial states denned by a = tc/A and a = n/20 respectively 
and an interacting case with k = 1.5, J = and a = 7r/20 shown in black. Both 
quantities, concurrence and purity, display oscillatory behaviour, with one frequency in 
the noninteracting case and two frequencies in the interacting case as can be verified in 
eqs. (fbil) and (fl5l) . One can also note that with interaction (black curve) the concurrence 
increases while the minimum value of purity is greater in contrast to the corresponding 
noninteracting case (blue-dashed curve). Similar behaviour in time domain has already 
been studied in other references like [21 El SI E] . 

The graphs in the time domain look pretty standard and this does not change if 
both interactions are present. It is therefore convenient to visualize the joint dynamics 
in a concurrence vs purity plane, the CP-Plane. Figure [2] a) shows the corresponding 
plane for the curves in figure Q] with the same colour code, but now the black curve is 
parametrized up to t = 20. In this plane we have plotted to guide the eye, a gray zone 
corresponding to the concurrence and purity combinations that can not be obtained 
in physical states and its lower frontier corresponds to the maximally entangled mixed 
states (MEMS), which for a given value of the purity maximize the concurrence [21J. 
The gray dashed line is defined by the Werner states pw = £z + (1 — 01 Bell) (Bell |, 
< f < 1 [211 [22]. 

One can note as well, that the dynamic of the interacting (black) case is enclosed 
by the noninteracting curves, the lower bound given by the blue curve with the same 
initial state as the black one, while the upper bound given by an initial state given by 
a = 7r/4. Perhaps the most important feature here, that one can not easily visualize 
in the time domain, is that for an initial bell state with a = tt/A and no interaction 
between the atoms, the curve (red) follows precisely, as we shall prove below, the one 
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that determines the mentioned MEMS. For this we need first to obtain the analytic 
solutions in the CP-plane. 

We take the explicit solutions in time in eqs. ( |T5l) . with n = J = 0, and invert 
them to find an explicit relation of the concurrence in terms of the purity. In this 
non-interacting case, concurrence is represented by up to two different curves in the 
CP-Plane 

Ci n) (P; a) = Max {o, |sin (2a) - f£\P)\ - n fi n) (P)} (16) 
with 7 n and f3 n as given in eq. ( TT4")) and with 



n 1 ± Jl + 7„ P- 1 

ft\p) = — v — 2^ 1 (17) 

In 

We find two separate cases: 

(i) For \ arcsin ( 3 J^3„ +1 ) < a < vr/4, the concurrence in the CP-plane is determined 
by the two curves: 

Cf(P;a), l-^<P< 7 " (1 + Sin(2Q))2 -sin(2.) 

In 4 

Ci n) (P;a), 1- — <P<1. (18) 

In 

(ii) Otherwise, the concurrence is determined only by the curve: 

C^(P;a), 7 " (1 + f (2a))2 -sin(2a)<P<l. (19) 

In figure [2] we show these solutions in the CP-Plane for n = and different values 
of a. The red curve shows the case when the starting state is the symmetric Bell state, 
a = tt/4. This solution has the explicit form 

CiV; V4) = ^ (l T v / 2P^I) (20) 

and it can be seen that in a certain region it coincides with the curve for the MEMS. In 
fact C^\P; tt/4:) coincides precisely with the curve of the MEMS for 5/9 < P < 1. The 
dashed blue curve represents the situation with an initial state determined by a pure 
but not fully entangled state. 

The behaviour for n > is qualitatively the same and when one takes as initial 
state the symmetric Bell state, a = 7r/4,the curves converges (n — > oo) to 



C M (P; tt/4) = Max {o, i (v / 24P^8 - l) } , i 



- < P < 1 
3 ~ ~ 

C<?\P; vr/4) = 0, \<P<\- (21) 

Actually in figure [3] we took n = 5 and the red curve is a very good approximation to 
C±°' ) . We note however that in the limit n — > oo this curve, which is actually an upper 
bound, lies below the werner curve. For finite n the curves Ci n) (P; 7r/4) intersects the 
werner curve in an additional point apart from C = 1. This means that there is a small 
region (hardly visible in figure [3]) above the werner curve that can be reached by the 
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Figure 2. (Colour online) CP— Plane for n — and an initially empty cavity. The 
red curve shows the behaviour in the CP-plane for two non-interacting atoms with an 
initial state determined by a — 7r/4 (See eq. l[7])) i.e. a maximally entangled pure 
state. The blue dashed curve shows the behaviour for non-interacting atoms with an 
initial pure, but not maximally entangled state a) a — tt/20 and b) a = n/10. In 
black, the curve for two interacting atoms with the same initial state as in the blue 
dashed curve and parametrized by time up to t = 20. a) k = 1.5 and J = 0. b) At = 1.5 
and J = 0.87. The gray area indicates CP combinations that can not be obtained 
in physical states and its lower frontier corresponds to the maximally entangled mixed 
states. The gray dashed line represents the Werner states. 




Figure 3. (Colour online) CP— Plane, same as fig. [2] but for n — 5. a) a = 7r/20 for 
the black and blue line and k = 5.7 and J — 0.2 for the black line, b) a = —ir/20 for 
the black and blue line and At = J = 5-^4 x 5 + 2 for the black line. 



dynamics. We do not write here explicit expressions for the dashed blue curves, as they 
can be obtained from equations ( TIB"]) and (jl9|) . 

Figures [2] and [3] show in black the curves for the interacting case for the same initial 
states as the blue dashed curves. One can note that in the CP-plane the curves now 
form a Lissajous-like figures with their frontier defined by the curves C±\P; a) for the 
starting value of a (lower frontier) and for a = tt/4 (upper frontier). Note that for 
increasing values of the difference of the interactions the curve in the CP-Plane does 
not fill the entire region enclosed by the curves Ci n) (P;«). The region filled by the 
black curve in figure |2] b) is smaller than in [2] a), because in [2] b) we use a larger value of 
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\k — J\. The minimum value of the purity can be calculated as -P m ; n = P(t = Tr/2u n ), 
this is the lower bound for P for the black curves in the figures [2)2) and(3]2). A separate 



the minimum value is P m i n = 1 — l/7 n , figure [2^). If both interactions have the same 
strength the curve will fill again the entire area, except in the case when there are 
commensurable frequencies in F and G, equation (fl4l) . That is the case of figure [3] b) 
where the black curve is closed. 

We have given closed solutions for the dynamics of two different TLA in a cavity 
interacting by dipole-dipole and Ising interaction. Many solvable models discussed for 
two TLA in a cavity belong to this wider class of exactly solvable models including 
a model for a Dirac Oscillator outside the realm of quantum optics [10J. The 
effectiveness of the general solution presented was displayed by calculating the evolution 
of concurrence and purity and fully determining the region of its evolution in a CP 
diagram in a particular, but interacting case. Interesting features appear when including 
both types of interactions. Intuitively one might think of less decoherence and a more 
robust entanglement with increasing interaction between the atoms. This is true if we 
have either of the interactions, but not necessarily if one takes interactions of similar 
strength. 

The parameter space of the model will be further explored in a full length paper. 
Interesting situations include placing one TLA outside the cavity or at a node of the 
mode and using different detuning as well as coherent or more complicated initial states. 

An interesting perspective would be to extend this technique to situations, where 
the couplings are chosen such that -rather than four dimensional spaces where exact 
solutions are available- we would have larger but finite spaces known in molecular 
physics as polyades, which are accessible to treatments with Lie algebraic [23j [24] and 
semi-classical [25], [26] techniques. 

Another worthwhile line of research may be to find an even more general class 
of solvable models including the one presented here and the ones using a pseudomode 
approach [3 [8] . 
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